Plane waves in a general Robertson- Walker background 
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We present an exact solution for the plane wave mode functions of a massless, minimally coupled 
scalar propagating in an arbitrary homogeneous, isotropic and spatially flat geometry. Our solution 
encompasses all previous solvable special cases such as de Sitter and power law expansion. Moreover, 
it can generate the mode functions for gravitons. We discuss some of the many applications that 
are now possible. 
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1. Introduction: The universe is homogeneous and 
isotropic on the largest observable scales . Recent pre- 
cision measurements of anisotropics in the cosmic mi- 
crowave background strongly support the idea that it is 
also spatially flat 0,0, E|. This is in any case the theoreti- 
cal prejudice of inflationary cosmology ||. The invariant 
element of such a universe can be simply expressed in 
co-moving coordinates, 



ds 2 



-dt 2 + a 2 (t)dx ■ dx 



(1) 



Although the scale factor a(t) is not measurable, its 
derivatives can be formed to give observables, 



H(t) 



«(*) 



1 -HH- 2 



(2) 



The Hubble parameter H(t) is the cosmological expan- 
sion rate while the deceleration parameter q(t) provides 
a measure of the rate at which this expansion is slowing 
down (q > —1) or speeding up (q < —1). 

It has long been realized that an expanding universe 
can result in particle creation 0. A particularly inter- 
esting and likely candidate for this effect is the massless, 
minimally coupled scalar, 



C 



(3) 



Because these particles are massless without classical 
conformal invariance, there can be an appreciable ampli- 
tude for virtual pairs to appear with wave lengths com- 
parable to H~ x . When this happens the pairs are pulled 
apart by the expansion of spacetime and there is physical 
particle creation 0. When a scalar potential is added 
to drive inflation, essentially this mechanism || is the 
leading candidate for explaining the primordial cosmo- 
logical perturbations which seeded structure formation 
and whose imprint is observed in the anisotropics of the 
cosmic microwave background [jj. 

An even more direct application derives from the ob- 
servation 0| that linearized gravitons obey the same 
equation of motion in the background (QJ. This leads to 



a fundamental tensor component in the primordial spec- 
trum of cosmological perturbations 0, 0, 0] . There 
is also the possibility that quantum gravitational back- 
reaction slows inflation . Although standard inflation- 
ary cosmology is realized via a scalar field, it is probably 
not natural to have interacting scalars which are much 
lighter than the Hubble parameter during inflation. In 
that respect the graviton can achieve all the advantages 
of inflation in a far more economical and natural way 

It is therefore frustrating that the scalar mode func- 
tions are only known for certain choices of a(t) 0|. The 
fact that typical inflationary backgrounds are not of this 
type means that either perturbative or numerical tech- 
niques must be used to pred ict the spectrum of cosmo- 
logical perturbations yjj, |l2|, . The absence of general 
mode solutions is most noticeable when back-reaction be- 
comes important. By working in de Sitter background 
one can see that the slowing effect eventually becomes 
strong [l4| , but it is not possible to reliably follow where 
the physics is trying to lead. 

The purpose of this paper is to solve the scalar mode 
equation for any scale factor a(t). Section 2 gives the 
equation and defines what we mean by the "mode func- 
tion" . In Section 3 we present a solution which is valid for 
any period during which the deceleration parameter © 
is constant. In Section 4 these solutions are extended to 
cover general evolution by representing it as a sequence 
of ever finer steps between regions of constant q(t). The 
transfer matrix which makes this possible is expressed as 
the path-ordered product of the exponential of a line in- 
tegral. In Section 5 we improve the constant deceleration 
solution by explicitly extracting the determinant of this 
matrix for general a(t). We also identify and normalize 
the positive frequency solution. Expansions for the ul- 
traviolet and infrared regimes are derived in Sections 6 
and 7, respectively. Our results are collected in Section 8 
and expressed in terms of standard parameters. We also 
discuss applications, including an improved result for the 
gravitational wave contribution to the power spectrum of 
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anisotropies in the cosmic microwave background. 

2. The scalar mode equation: We take advantage of 
spatial translation invariance to Fourier transform in the 
(D— l)-dimensional spatial coordinates, 



<p(t, k) = / dr~ht e^Xf, x) 



The equation of motion of (p(t, k) is, 
ip+ (D-l)Hp+ k 2 a' 2 lp 







(4) 



(5) 



The middle term, (D—l)Hlp, is responsible for the phe- 
nomenon of "Hubble friction" in which the expansion of 
spacetime retards the scalar's evolution. 

Our work was motivated by a recent paper of Finelli, 
Marozzi, Vacca and Venturi who studied the case 
of a minimally coupled scalar with arbitrary mass. We 
follow them in the standard step of extracting a factor of 
a ~ to eliminate the Hubble friction, 

<p(t,k) = a^^u(t,k)a(k)+u*(t,k)a t (-k)^ . (6) 

The quantities a(k) and cr(— k) are time independent 
quantum operators. What remains is the C-number 
mode function u(t, k) and its conjugate. It obeys, 

u(t, k) + H 2 [(^) 2 - - {^) 2 } u(t, k) = . (7) 

To fix the normalization, suppose the creation and an- 
nihilation operators commute canonically, 



a(fc),a t (fc') = ^^S^ik-k') 



(8) 



Taking the time derivative of © gives, 

D-l 



2tt 



S^ik + k') 



x(u(t,k)u*(t,k) - u*(t,k)u(t,k)\ . (9) 

Since the momentum canonically conjugate to ip(t,x) is 
a- D_1 (<)(yj(t, x), we must also have, 



2?r 



D-l 



ftt,k)Mt,i?) =l i^itj) ^{k + k'). (io) 



Comparing JSJ with l|l(J|) determines the Wronskian of 
the mode function, 



u(t, k) u*{t, k) — u*(t, k) u(t, k) = i 



(11) 



3. The solution for constant deceleration: Just as in 
|l(j|. we organize the dependence upon time and wave 
number into two dimcnsionlcss parameters, 



x(t,k) 



k 
~a~H 



y(t) 



H 

IP 



(12) 



The meaning of x(t, k) is the physical wave number in 
Hubble units. Super- horizon modes have x < 1 whereas 
sub- horizon modes have x > 1. The quantity y(t) is 
related to the deceleration parameter, q(t) = — 1 + y{t). 
The Weak Energy Condition implies y(t) > 0. Slow roll 
inflation is characterized by y{t) -C 1. 

Note that y(t) is constant whenever a(t) is a power 
of any linear function of t. For example, a radiation 
dominated universe has y — 2; the result for a matter 
dominated universe is y — |; and curvature domination 
gives y = 1. Of course J7J is trivially reducible to Bessel's 
equation whenever the scale factor obeys a power law, but 
the order of the necessary Bessel functions varies with the 
power. This is why it has been so difficult to find a gen- 
eral solution. In [lfj, however, it was recognized that the 
general solution should involve a Bessel function whose 
order depends upon y(t). They worked only to first order 
in y but it is straightforward to extend their result to all 
orders. The required generalization is, 



H(x,y) = x^H^(Xx) 



(13) 



where the parameters are the following functions of y, 

^=T^r. > v =—, — — > x =t—.- (14) 



1 



y 



l 



y 



Although the parameter x is always positive, A passes 
from +oo to — oo as y evolves from 1~ to 1 + , which sig- 
nals the end of inflation. This does not lead to wild 
oscillations in 7i(x,y) because v — \ + (y — 1)A makes 
the same transition, and the factor of 1/T(n + v + 1) in 
the series expansion of the Hankel function compensates 
the growth of A. However, the possibility of A becom- 
ing negative does raise the issue of how we define (Ax) 1 '. 
Throughout this paper we shall understand negative A to 
mean e+"|A|, so that (-\\\x) v = {\\\x) v . 

To see that (|13|) solves (J7J) for constant y, first compute 
the (general) time derivatives of x(t, k) and y(t), 



H 



H 



k 

a 


] + W 2 _ 


= -Hx(l - 


Hk 


H 


H 2 H 


a 


+ 2 

H 4 H 3 




= H 2 x 





(15) 



(16) 



(17) 



Now use Bessel's equation to derive an identity for the 
second x derivative of Ti(x, y), 



.XX 



2/1-1 

x 

7~L,x 
x(l - y) 



n.. 



D-l 

2 



V 2 - /i 2 
X 2 



x 2 (l — y) 2 



n , (is) 



H . (19) 
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When y = we have, 



H 

n 



y=o 



xH. 



— x T~~t x H~ 7~t,xx ) 



(20) 
(21) 



= H 2 x{l-y)H, x + H 2 x 2 {l-y) 2 H, xx , (22) 
= -H 2 [x 2 + S^y- (%i) 2 ] H • (23) 

4- The transfer matrix: Let us follow Grishchuk [l7j ] 
in joining solutions between two periods of constant y. 
Suppose there is a sudden transition at t = tx, that y(t) = 
yi for t < tx, and that y(t) = y<% for t>t\. The solution 
during each period must be a linear combination of the 
appropriate constant y solutions, 

ux{t,k) = cxH{x,yx) + dxH*{x, yi ) Vt<t x , (24) 
u 2 (t,k) = C2H(x,y2)+d2H*(x,y 2 ) Vt>i x . (25) 

Suppose further that the scale factor is continuous at the 
transition. Since the differential equation has no factors 
of S'(t — tx) its solution must also be continuous at tx, 



C2^12 + d 2 Hl 2 = c{H\\ + dxH*i , 



(26) 



where we have compressed the notation by making the 
definition, Tiij = Ti(x(ti), yj). Suppose finally that the 
Hubble constant is continuous. Since the differential 
equation has no factors of 8(t — tx) the first derivative 
of its solution must be continuous at the transition, 



(1-2/2) 



c-i7~L\i,x + d 2 Tt 12x 

= (i-yi) 



cxHxx + dxH* n<x . (27) 



We represent the linear transformation between the 
"before" and "after" combination coefficients as a 2 x 2 
matrix T(ii), 



= T(ti 



(28) 



Solving for its components is a straightforward exercise 
in linear algebra, 



T 11 ^) = [T 22 (tx) 
rtxxrL l2 



1 

w 2 



T 21 (h) = [T 12 (tx) 
1 

~ W 2 



— T~txxT~ix2,x 



1-2/1 

1-2/2 
1-2/2 



riXX,xrtx2 



T~ixx,xHx2 



(29) 



(30) 



Here W2 = W\x(tx), 2/2 J and y) is the Wronskian 

of our constant y solution l|13|) , 



W{x,y)=HH* x -H <x H* 



4? 



a; 2 ^~ 1 sgn(A) . (31) 



Now consider a sequence of such transitions at times 
tx < £2 < £3 < • ■ ■ The solution for t > t n _i can be 
expressed as the inner product, through a product of 
transformation matrices T(tfc), between the row vector 
(H , 7i*) and the column vector composed of the origi- 
nal combination coefficients, 



u n (t,k) = (H(x,y n ) , H*{x,y n )} [] T (*n-fe) ( j 

fc=i ^ 1 

(32) 

We can approach continuous evolution in y{t) by repre- 
senting any fixed interval as a series of more and more 
transitions of ever smaller magnitude. In the infinitesi- 
mal limit the product of the T matrices exponentiates. 

To demonstrate exponentiation consider the transition 
at t = tx for yx = y{tx) ~ §Ay and y 2 = y(h) + \Ay. 
The Wronskian and the constant y solutions can all be 
expanded in powers of Ay, 



(33) 
(34) 
(35) 
(36) 



" ) 


= 1 + 


Ay 




l-2/ + ' 


2/2/ 




1 


1 


W, y Ay 




W 


W 2 2 


Tin 


= H - 




7~ix2 


= H + 


H -y— 4 



Expanding the coefficients of T(ii) to first order gives, 



T n (ix) = l+c(a;(ti),|/(t 1 ))Ay + 
T 21 (h) = + t(x(tx),y(tx))Ay + . 

where we define, 



— Tt,yTt* x + TC^yTi* 

~H,yH,x — 7~L. xy Tl + 



rt, x l<- 

~ 1-2/ 
'H.x'H 

1-2/. 



(37) 
(38) 

,(39) 
(40) 



Taking the infinitesimal limit gives the path-ordered 
product of the exponential of the line integral of the fol- 
lowing matrix: 



A(t,k)=y(t) 



C(x,y) C(x,y) 



J(x,y) C(x,y) 
We define the transfer matrix as, 

M(t,U,k) = pjexp J dt'A{t',k) 



(41) 



(42) 



dt, 



xA(tx,k)---A(t n ,k). (43) 
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Here is the initial time. The full solution assumes the 
form, 

u(t,k)=(n(x,y) , H*(x,y))M(t,U,k)(° \ . (44) 

Explicitly verifying that (|44|l obeys Q is facilitated by 
two identities implied by the Wronskian, 



(H + £H* = -H 



0~L,x + iH* 



T~L,x 



(45) 
(46) 



Another very useful fact is that the time derivative of the 
transfer matrix is proportional to it, 



M(i, t h k) = A(t, k) M(t, U, k) 



(47) 



Finally, we recall that time derivatives of TL (x(t, k), y(t) 
follow from the chain rule, 



H(x, y) = xTL x + yH, % 



(48) 



With these facts in hand it is straightforward to com- 
pute the first derivative, 

u(t,k) = [(ft , n*^j + (n , w*)a]mQ, (49) 

d) ■ ^ 

The second derivative requires only a little more effort, 



= x(H 



H*)M 



u(t,k)= x[H lX , H* x ) + x[H lX , H* x 



M 



xH, x +x*H, xx + - — - , c.c. M . 
■H 2 (x 2 + ^y-(^) 2 )u(t,k) . 



(51) 

(52) 
(53) 



In the last step we have substituted expressions (|16I17|) 
for x and x, and we have exploited relation l|19|) for 7i. xx . 

5. Evaluating the transfer matrix: During most peri- 
ods of interesting cosmologies one has \y/H\ <C 1. Since 
M vanishes with y it is reasonable to expect that the 
mode functions derive most of their time dependence 
from the constant deceleration solutions, Tt(x,y) and 
TL*{x,y). This is one reason our solution is effective in 
spite of the cumbersome definition i|39l43[l of the transfer 
matrix. The role of M(t, U,k) is to slowly mix positive 
and negative frequency solutions as cosmological particle 
production progresses. It must be a nonlocal summation 
over the past. 

The nonlocal character of M(t, U, k) does not preclude 
learning a great deal about it. We begin by recalling 



that H(x,y) = x tJ, Hi 1 \\x) with fi — ^ — and v = 
h + (y — 1)A. Now apply the chain rule to show, 

H, x = pix' J - 1 Hi l \Xx) + Xx tt H l ^ l {\x) , (54) 
H, y = fi'\n(x)x^Hl 1 \Xx) 

+X'x» +1 Hi 1) '(Xx) + v'x^d„Hl 1 \Xx) . (55) 

Here a prime means differentiation with repsect to the 
function's argument. Upon substitution in the definition 
(|39|1 for y) and making some tedious rearrangements 
we find, 



C= l -XHn{x)+ l -X 



ITV 

T 



\X\{^El(Xx)+F u (Xx) + (v-i)G v (Xx)}, (56) 



where E v (z), F u (z) and G v {z) are the following real prod- 
ucts of Hankcl functions, 



E v {z) 
FAz) 
G v {z) 



H?\z) 



dz HP(z)\\+(l-^)\\HP(z) 



(57) 



, (58) 



-d v d.H, 



f){z)(H^{z)] 



(59) 



It might be thought that these three functions require 
separate expansions but the last two actually follow from 
E v (z) through the identities, 



F' v {z) = 2E v {z) 



G' v {z) = ~E v (z) . (60) 



The same procedure can be used to bring £(a;, y) to a 
closely related form, 

i = f|A|{~4(Az) + F„{Xx) + §)&,(»}. 

(61) 

In this case each function is the product of the same two 
Hankel functions, 



u {z) = z(H£\z))\ 



r v {z) 

GAz) 



(d z HU(z))+(l-^)(HP(z)y 
d u Hl 1 \z)d z Hl 1 \z) 

-d v d z H^{z)H^{z) 



(62) 



, (63) 



(64) 
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These functions are not purely real, unlike those for £. 
However, they still bear the same relation to one another, 



T' v {z) = 2£ v {z) 



G' v {z) = —Euiz) . (65) 

z A 



It is straightforward to obtain explicit series expan- 
sions for C,(x,y) and £(x,y). One begins by expressing 
E v (z) and £ v (z) in terms of Bessel functions, 

sm (vir) i. 

-2cos{v<k)J- v {z)J v {z) + zJt{z)\ , (66) 
E ^ z ) = —in — A- zJ -A z ) 

sm {vk) l - 

+2e-* w J_ v («)J v («) - e- 2ira z^(z)} .(67) 
Now exploit relations 1)60(1 and (|65[) to show, 



4 sin (W) 



[b v {\X\x)- 



+4sin 2 (^7r)6'(-A)+cos(^7r)c l ,(|A|a;) + d t/ (|A|a;)|, (68) 
~i-k{v — \)\X\~ 



4 sin (j/7t) 



{6„(Ax) - 2 
'"c„(Ax)+e- 2 ^d,(A 2 :)}, (69) 



where the various coefficient functions are 

K{z) 



i » (-i^r^-iz-D^^n-!/)- 1 (?o) 



d v {z)- 



2^pa 4^ r(n)r(n - 1/ + i)r(n — 2^ + 1 

4 . r 
sin(^7r) — 1 — m(|z) 

1 ~ (_l)«r(n- i)z 2n n-! 
~ ^ r(n+i/)r(n + l)r(n-i/+l)' 



(71) 



1 ~ (-l)»r(n+ i /-i)z 2 "+ 2 -(n + ^)- 1 , 

w„ , o I TvFTI rTT ' ( 72 ) 



20F ^ r(n + 2u)V(n + v + l)T(n + 1) 



Recall that ip(z) = T'(z)/T(z). 

It is illuminating to express the matrix A(i, fc) (|41|) in 
terms of the real and imaginary parts of £ and £, 

A(t, fc) =yRe(C)I+?/(Re(£)cr 1 +Im(Ocr 2 +iIm(C)cr 3 ) . (73) 



i-^-ln(A.T A a)l + A(t,k) 



(74) 



The unit matrix I commutes with all other matrices so we 
can factor this part out of path-ordered product. Since 
its coefficient is a total derivative we can also give an 
explicit form for the contribution it makes to M(t, ti, k), 



M(t,U,k) 



Xx x a 



A M(t,t h k). (75) 



The subscript "i" in this and subsequent expressions 
refers to quantities evaluated at the initial time t = ti, 
for example, = a(tj). It is worth emphasizing that 
extracting the determinant of the transfer matrix already 
gives an improvement on the previous constant y solu- 
tions. 

The residual transfer matrix M(t, ti, fc), 



M(t,t h k) = P< exp 



dt'A(t',k) 



(76) 



has unit determinant and can be represented in terms 
of a nonnegative real number p(t,ti,k) and the phases 
r(t, U, fc) and u>(t, U, fc), 



M(t,U,k) = 



2„rr 



pe 



pe" 



V 1 + P 2e 



(77) 



We can use these parameters to express the fully normal- 
ized positive frequency solution. First note that most of 
the time dependence of the full solution resides in 7t(x, y) 
and Ti*(x,y). Next note that the argument of the Han- 
kel function in TL{x,y) is Ax. From the time derivative 
of this quantity, 



) = -Hx 






H 



(78) 



we see that Xx is typically a decreasing function of time. 
It is therefore reasonable to choose the positive frequency 
solution to be initially pure H(x,y). With our normal- 
ization this means that the full solution is, 



u(t,k) = sj { y/TT?e iT H(x, y) + pe^H* (x, y) } . 

(79) 

6. The ultraviolet regime: One can obtain a simple ap- 
proximate form in the ultraviolet regime of x 1. Re- 
call that the part of £(x, y) which appears in the residual 
transfer matrix is, 

l{x,y)=~x{^E' v {Xx) +F v {Xx) + {v- ^)G v {Xx)]. 

(80) 

Recall also the well-known asymptotic expansion for the 
Hankel function of the first kind, 

\i Try — ' nil 1 7/_n X± \ ZlZ ' 



(82) 



With this to set the integration constant we can use re- 
lations (|60|l to derive the following exact results, 



\e' v {z)+F v {z) = -z-2{v 2 -^)z 

l 7T 



Iz'^fi, (83) 



G v (z) = -2 + 2u / dz 



,E v (z') 



(84) 
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Now substitute l|81[l and integrate termwise to obtain an 
asymptotic expansion for ((x,y), 

C{x, y) = -iX 2 x + i$(u- ±)A 

i{v-\)X ~ r(n+l)r(^+n + l) (Ax)' 2 "- 1 
2 ^ h {n+\)\T{y-n-\) (2n + 1) ' 1 j 

It is simplest to express the asymptotic expansion for 
£(x,y) in terms of the function h v (z) defined in 



+«[/i;(2)] 2 +ft,^(«)/»u«)-W(«)^(«)} • (86) 

From the expansion of h v (z) we see that £(x, y) is neg- 
ligible in the ultraviolet, 



o((Ax)- 3 )}. 

(87) 



Ignoring £(x, y) defines the ultraviolet regime, 

1 



A uv (f, fe) = y((x,y) 



-1 



(88) 



In the notation of (|79|l the ultraviolet regime is charac- 
terized by, 



p uv (t,U,k) = , 
T uv (t,ti,k) = J dt'y((x,y) 



(89) 
(90) 



Since p uv = the value of lu uv is irrelevant. 

In the extreme ultraviolet limit we need only keep the 
first two terms from the expansion of £(x, y), 



7"uv (t: ti > ^) 



dt'y -\ 2 x + § (v- ±)A 



(91) 



-\x + f(v-i) 



dt: 



The integral proportional to — k on the last line is just the 
conformal time interval r\ — rji. Substituting as well the 
asymptotic expansion of the Hankel function in 7i(x,y) 
gives the following leading ultraviolet result, 



u m (t,k) — > \/ — e 



fl W „— ifc(??-J?i)+ix» 



(93) 



where x(^i^) = A(t) a;(t, fc)— [^(i) + |]f ■ Recovering the 
correct correspondence limit is an important check on the 
consistency of the formalism. 

7. The infrared regime: For phenomenology it is the 
infrared regime (i<1) that is more useful. In this case 
it is best to decompose the Hankel functions into Bessel 
functions of positive and negative order, 



H(x,y) = icsc(^7r)a: M 



e- w *J„{\x) - J_(\x) . (94) 



The conjugate solution takes slightly different forms for A 
positive (inflation) and negative (subluminal expansion), 



Ti.*(x, y) = icsc{vir)x^ 



= ii/7rsgn(A) 



J v {Xx) 



+e 2i " n0( - x) J^(Xx)\. (95) 



We seek a matrix 0(t) that effects the change of basis 
from H and H* to solutions based upon J_„ and J„, 



-j=(H(x,y),H*(x,y))o(t) 



x M ^— icsc(^7r) J_„ (Xx) , J u {Xx) 



(96) 



For positive A (inflation) the required matrix and its in- 
verse are, 



, . _ J_ / l-icot(i/7r) 1 
W ~ V2 ^-l-icot(w) 1 



O 1 = 



w i 

V l+«cot(^7r) 1 — icot(z/7r) 



For A negative (subluminal expansion) we have, 



, J_ ( l+«cot(^7r) 

W ~ V2 V,l+*C0t(^7T) 



-e 
-1 



O 1 



1 / 1 e 2il/7r 
\ l+icot(i/7r) — 1 — icot(^7r) 



(97) 
(98) 

(99) 
.(100) 



In the new basis this matrix is, 

Mih^k) = 0-^ 2 )M(t 2 , < l7 *)0(ti, fc) , (101) 



= P < exp 



/c) 



(102) 



The new exponent matrix is related to the old as follows, 

A{t, k) = Or\t)0{t) + 0~\t)A(t, k)0(t) . (103) 

This is just an exercise in matrix multiplication. For 
either sign of A the result is, 



^ _ 7T ^ ( csc(vTt)c v {\x) 



—2id v {Xx) 
4 V— 2zcsc 2 (z/7r)&j,(Aiz;) — csc(i/7r)c„(Ax 

j(t,k) -iS(t,k) 



-tf3(t,k) - 7 (t,fc) 
For small x the coefficients of A go like 

(3 ~ x 2 - 2u , 7 - ln(x) , J 



2,i 



(104) 
(105) 

(106) 



Although [3(t, k) makes the dominant contribution to the 
residual transfer matrix for small x, it does not make 
the dominant contribution to the mode function. This is 
because the only way a factor of f3 can multiply the J- v , 
rather than the J„, is if it accompanies a factor of 5. To 
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see this, first consider the basis functions times a single 
factor of A, 

icsc(z^7r) 3- v , jJjA 
= (— icsc(W)7J_„— if3J v ,— csc(^7r)<5 J_„— 7 J^j . (107) 

The small x behavior of 7J_„ is ~ \n(x)x~ l/ , whereas 
(3J U ~ x 2 ~" is down by a; 2 . Now consider two factors 
of A for any times during which x <C 1, 



^4(t2, A)^4(ti, ft): 



l2li~faPi -i[72<5i-<5 2 7i] 
-«[/?27i-72/3i] 727i-/?2<5i 



(108) 

Since /3<5 ~ x 2 , these terms are negligible with respect to 
7 2 ~ ln 2 (x). 

To the same proportional error of order x 2 C 1 we can 
easily obtain an explicit expression for the n-fold product, 



A(t n ,k)---A(tx,k) - 

lnln-l ' ' ' 71 

-«[/3„7„_i •••71- 



-i[ln ■ ■ -72<5l-- ■ ■] 
(-l)™7r ! 7r ! -l-"7l 

In this expression the off diagonal elements are, 



(109) 



-i[ln ■ ■ '12 fa 



In ■ ■ 'llfa 



-In ' ■ -73<527l + - ■ --(-l)"<5n7n-l 1 ' ' H ) (HO) 



-i[/3„7„_i • • -7i-. . 



0nln-l ' ' ' 71 



-^fnSn^-lln-2 ' ' ' 7l + ■ ■ •-(-l)™7n ' "1201 • (HI) 

This form can be exponentiated, 

Mr(t2,£i,ft) — 
where the various components are, 



r(* 2 ,*i,ft) -iA(t 2 ,h,k) 
-iB(t 2 ,h,k) r-\t 2 ,h,k) 



(112) 



r(*2, h, k) = exp 



dtj(t, k) 



(113) 



The leading term in (|117|l is a total derivative, 



i(t,k) = -m 



a\xy 



+ 0(x 2 ), (119) 



so we can obtain explicit expressions for T, B and A, up 
to correction factors of order x 2 , 



T(t2,ti,k) 
B(ta,ti,fc) 



h 2 \^2X2Y 2 - 1 i 



i> 2 ) 
i> 2 ) 



< 1 (|Ai^i) l/1 - 

i>i) 



(120) 



a 2 2 \\\2X2Y^a? \\Mx x y^ 



dt- 



tXQ 



D-2 



A(t 2 ,ti,fc) 



2n(D-2)J ti (v-lf 
af-\\\ 2 x 2 y^ af^ilXxXxY-- 



(121) 



i> 2 ) 



2n(D-2) 



dt 



v x 



(122) 



It is useful to compare the elements of the same column 
of the infrared transfer matrix (|112|l . The ratio of the 21 
and 11 components is, 



-iB(t 2 ,h,k) -iT 2 {v 2 ) (I 



r(* 2 ,*i,*) 



2n(D-2) 



x 2 



{^2X 2 



B(t 2 ,t u k) = / dtr-\t 2 ,t,k)l3(t,k)r(t,ti,k), (114) 

A{t 2 ,h,k) = ^ dtT{t 2 ,t,k)5(t,k)T \t,h,k). (115) -iA(t 2 ,t u k) -iirx^X^) 2 ^- 1 [*> v H 2 fa 2 ^-i 



i> H 2 ( a \ D ~ 3 / 

/'FIF^U - (123) 



In the infrared regime we know £2 *C 1. The expansion 
of the universe implies a(t)/a 2 < 1, and the weak energy 
condition tells us H 2 /H(t) < 1. During inflation the 
parameter v{t) is positive (y > -^p") so the factor a; 2 , -2 " 2 
can win out over the other terms and make —iB(t 2 , t\, k) 
dominate T(t 2 ,t\,k). After inflation the parameter v(t) 
changes sign and —iB(t 2 , t\, k) is insignificant compared 
with T(t 2 ,ti,k). (The only exception is when v 2 passes 
through a negative integer, at which point T(v 2 ) diverges 
and the infrared basis becomes degenerate.) The ratio of 
the 12 and 22 components is, 



Because <|112[) is only valid up to factors of order x 2 it 
is superfluous to evaluate the components using anything 
but the leading small x forms for (3, 7 and 5, 

j(t, k) = z>[hi(-Ax) + 1 - VM] + 0{x 2 ) , (117) 
8(t, k) = 1 \^ X? " + 0{x 2 ^ 2 ) . (118) 



/02 

V a 



■ 



r~\t 2 ,ti,k) 2{D-2)T 2 ( V2 ) J t ;"v 2 H\a) 

(124) 

Some of these factors are less than one, others are greater, 
and the magnitude of the ratio is not clear. 

Recall the series expansion of the Bessel function, 



_ ^ { _ 1)n{ l z) 2n +V 

u(Z, ~ ^ n\T(n + l + u) 

71—0 K 7 



(125) 



2i/-|r 2 (i/-i) 



From this we see that the linearly independent solutions 
appropriate to the far infrared consist of a "big" mode 
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function u (t,t±,k), based on J_„(Xx), and a "small" 
mode function u + (t,ti, k), based on J^(Xx), 



u (t,tl,k) 



u + (t, ti, k) 



ira Xx ( 

VttI' 



icsc(;/7r) J- v (Xx) , 



J v (Xx))M(t,ti,k) 



na Xx 
~k ~2 



(—icsc(vTr)J- u (Xx) 



J v (Xx))M(t,ti,k) 



(126) 



(127) 



It is straightforward to verify that these solutions obey, 

vT (t, ti, k)ii + (t, ti, k) — ii~(t, t\, k)u + (t, t\, k) = i . 

(128) 

If x <C 1 between t% and t, we need only keep the n = 
terms of the Bessel functions (I125|l . We can also employ 
the infrared simplifications (|112I115J) of the transfer ma- 
trix. Even though B(t, ti, k) can be larger than T(t, t%,k), 
its contribution to u~(t, ti, k) is still insignificant because 
it multiplies the J„, 



u (t,ti,k)- 



na Xx 
~k ~~2 



D-l 

-ia 2 



csc(i/7r) J-„(Xx)r(t, ti, k) 
+J u (Xx)B(t,ti,k) , (129) 
I>i) 



'irk 



\\x lXl y^ 



(130) 



Most of our limiting form H130|) consists of constants; its 
time dependence derives entirely from the factor of a~ ~ . 
This corresponds to a constant scalar field, and we see 
from JHJ) that tp(t, k) = constant is indeed a solution for 
the infrared limit of k = 0. 

Even in the infrared limit both terms contribute to 
u + (t,ti,k), 



u + (t, ti,k)- 



na Xx 
T T 



, — p-i a 
irk a 2 — 



-csc(^7r) J-v(Xx)A(t 7 ti,k) 

+J u (Xx)T-\t,t 1 ,k)\,(131) 



1 



2D-4 



dt' 



X 



v 2 Ha D-i 2vHa D ~ l 



(132) 



It is straightforward to partially integrate once, 



dt' 



1 



v 2 Ha 



D-l 



vHa 



D-l 



dt' 



va 



D-l 



From the identities, 

D-i-y 

v 
1 



2 — iQ 

x 



2 = (D — 2) 
v v 



(133) 

(134) 
(135) 



we see that a second partial integration brings the inte- 
gral to the form, 



■1 



2D-4 



dt' 



-X 



dt' 



(136) 

Substitution in (|132fl gives the infrared limit for the small 
mode function, 



f— a} \\X x xxf^ 
'irk a 2 - 



I>i) 



Ai 



2v 1 H 1 a°- 1 



dt' 
7d~i 



(137) 



Note that the infrared limits I|13UI137[) of our solutions 
retain their canonical normalization ^ 128(1 . 

The two independent solutions of the scalar field equa- 
tion (jsj) in the infrared limit of k — are, 



lp(t,0) 



1 



ip(t,6) 



dt' 



(138) 



We have seen that u~ (t,ti,k) becomes proportional to 
the constant solution, but the second solution in 1|138|) 
also contains a constant term from its dependence upon 
the lower limit of integration. Since the integrand falls as 
the universe expands, care must be taken to isolate the 
small, time-dependent term from the potentially much 
larger constant. This is something our formalism does 
automatically. For if y(t) is constant for t > t\ we can 
solve for the Hubble parameter and the scale factor, 



H{t) 
a(t) 



Hi 



1 + yxH^t-ti) 



0.1 



l + Vifri(i-ii) 



(139) 
(140) 



In this case the integral can be explicitly evaluated, 



dt' 



■.D-l 



Ai 



Ai 



y=vi 



2vxHa 



D-l 



2v\H\a x 



D-l 



(141) 



Hence the constant on the second line of our infrared 
form i|137|) for u + (t, ti, k) cancels the large constant from 
the lower limit, up to very small terms coming from the 
weak time dependence of y(t). So our second solution 
is indeed dominated by the small, time-dependent, lower 
limit term. 

It remains only to express the original mode solution 
(|79|l in the infrared basis, 



i(t,k) = -^=(u (t,ti,k) , u + (t,ti, 



k) 



xM(ti,U,k) 



l+icot(^i7r) 



(142) 



It is reasonable for ti to be the time of first horizon cross- 
ing, that is, x(ti,k) — 1. Even though x(t, k) is just 
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becoming small at this time, the expansion is so rapid 
during inflation that we have x <C 1 almost immediately 
afterwards. Because the transfer matrix up to this time 
is of order one, only u~(t, t±, k) contributes significantly, 



u iT (t,k)- 



-ia— Y(\-v)J_ v {\x) 
(IA*)-" 



2k 



(-^) 2 Ci(k)d(k). (143) 



We have expressed itj r with the standard normaliza- 
tion times two correction factors. The correction that 
depends upon the system at first horizon crossing is, 



Ci(fc)= f \ 



(144) 



The other correction factor depends upon evolution up 
to this point, 



C i {k) = M 11 (t 1 ,t i ,k) + 



-M 12 (h,ti,k). (145) 



sin(^7r) 

Further corrections to l| 143(1 are down by a factor 
x 2 (t',k) for some t' between first and second horizon 
crossings. (Second horizon crossing is the time ti 1 after 
inflation, at which x(t2,k) — 1.) Before second horizon 
crossing we need retain only the leading term in the se- 
ries expansion of the Bessel function ((125(1 . but the form 
((143(1 should remain valid to arbitrarily late times pro- 
vided that y(t) does not experience violent evolution at 
late times. This is a very safe assumption given what we 
know about late time cosmology. 

8. Summary and discussion: This has been a long and 
technical series of derivations, and it was not clear at 
the beginning what would be the most economical form 
in which to express the final answer. Now that we are 
done, it is well to summarize the important results in a 
simple way. It is also desirable to replace the parameter 
y{t) = —H/H 2 , and the quantities A and // which are 
derived from it, with the familiar deceleration parameter, 
q(t) = —1 + y{t). The only exception we make is the 
parameter v = i — (y — 1)^, which appears in too many 
indices to be conveniently expunged. 

Our mode solution is specified in terms of parameters 
constructed from the scale factor and its derivatives, 



x(t, k) 



k 
~a~H 



H 



and v{t) 



12-1 

2 x 



(146) 

The solution to equation JJJ) which is positive frequency 
at the initial time t = U takes the form of a row vector 
of Bessel functions multiplied into a transfer matrix, 



u (^ k )=\r^:\-^z) 2 (- icsc ( I/7r )- 7 -^(-f)) J ^(-f) 



2k V 2q 



X ^^^ fc )(l + l cotW,h])- (147) 



The transfer matrix Ai(t, ti, k) is the time-ordered prod- 
uct of the exponential of a line integral, 



M(t,U,k) = exp 



dt'A(t',k) 



(148) 



dh dt 2 . 



dt n A(h,k) ■■■ A(t n ,k) (149) 



The exponent matrix A(t, k) vanishes whenever q(t) is 
constant. It has the form, 



A(t, k)=-v 



_7T / csc(^ 7 r)c I/ (-f) 



-2*d„(-f) 



-2icsc 2 (i/7r) &„(-§) -csc(vk)c v {-z)J ' 



(150) 



where the various coefficient functions are, 



^ 1 ^ (_l ) n r(n _ 1/ _l )z 2^ (n _ y) -l ^ c ^ 

vW ~2^f£j T{n)T{n-v+l)T{n-2v+l) { ° ' 

4 r -I 

c v {z) = sin(j/7r) ip{v) — 1 — ln(iz) 

7T L 

~t r ( n+I/ ) r ( n + 1 ) r ( n - I/ + 1 ) 

1 g (-l)"r(n+^-l)z 2 "+ 2 -(n + ^)- 1 (i53) 



d v {z)- 



2^ 



n=0 



T{n + 2v)T{n+v + l)T{n+l) 



In the ultraviolet (x(t, k) 3> 1) our general solution 
((147(1 approaches, 



o(t) 



(154) 



werex(i,^) = A(t)a;(*,fc)-[i/(t)+§]§. This is just the flat 
space solution with a trivial redshift factor. Corrections 
to are of order l/x times ((154(1 . 

It is common to assume all modes are originally sub- 
horizon (x(ti,k) > 1). During inflation x(t, k) falls, typ- 
ically exponentially. This carries more and more modes 
into the infrared regime of x(t, k) <C 1. First horizon 
crossing is the time, t±(k), at which x(ti,k) — 1. After 
this time our general solution 11471) takes the form, 



r (t,fc)- 



r(i-^)j_ !/ (-f) 



2/v 



__x_\—v 
2q> 



x(^i) 2 'c^dik). (155) 



Corrections are of order x 2 {t' , k) times this result, where 
t' is some time for which x(t' , k) -C 1. 

The factor C\(k) depends upon the state of the system 
at first horizon crossing, 



(156) 
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The factor Cj(fe) depends upon evolution up to this time, 



C i (k) = M 11 (t 1 ,t l ,k) + 



sin(^7r) 



M 12 (tut h k). (157) 



Since both y(t) and y(t) are small during inflation, we 
can obtain an excellent approximation for Cj(fe) by simply 
expanding the transfer matrix, 



M lx {t u ti,k) = 1 + / dtj(t,k) 



ti 



dt I dt' 



7(i,fc)7(i',fc)-<5(t,fe)/3(t',fe) 



M (ti,ti,k) = i y dtS{t,k) 



-ij dtj dt' 



j(t,k)5(t',k) - 6{t,k)j{t',k) 



(158) 



(159) 



For completeness we remind the reader of the relations, 



0{t,k) 
l{t,k) 



2 sin (utt) 

717> 



4sin(^7r' 



Cv(-f), 



<5(t,fc) = — 



(160) 
(161) 
(162) 



The formalism we have just summarized has many ap- 
plications. One of these is to derive an improved estimate 
for the graviton contribution to the power spectrum of 
anisotropies in the cosmic microwave background. We 
shall present the derivation elsewhere but the result 
is, 



V h {k) = GHi{k)\\Ci{k)\VCi{k) . 



(163) 



The standard result is GHf, so the factors Ci(fc) and 
C\(k) represent improvements. Because the literature 
abounds with different conventions for this quantity we 
correspond Vh{k) below to the symbol 5h(k) used by 
Mukhanov, Feldman and Brandenberger |ll|. to the vari- 
able V g {k) used by Liddle and Lv th 1121 . and to the quan- 
tity Ay(fc) used by Lidsey et al. 13], 

Q77- nr ZOTT 

V h (k) = -\\5 h (kW = YQ P g (k) = — Al(k) . (164) 

The most unambiguous definition of Vh(k) is given by the 
(purely gravitational wave contribution to the) correla- 
tion function between temperature fluctuations observed 
from directions ef and €2, 



ATfl(ei) AT R {e 2 ) 



Tr T r 

xe*(e 2 ,k)ele(e^ 



1 



(165) 



In this formula Hy = Sij — kikj is the transverse projec- 
tion matrix coming from the sum over graviton polariza- 
tion tensors. The angular factor is, 



6(e, jfe) 



a 2i9 x 



47T 



2-W z 



1- 



ln 



1- 



1- 



!)H}' 



(166) 

where xq = x(to, k) is the physical wave number in cur- 
rent Hubble units and we define 6 = k ■ e. 

There are many other applications. One particularly 
interesting question we now have the technology to an- 
swer is the response to multiple periods of inflation. 
Whereas the weak energy condition implies that the Hub- 
ble parameter can only decrease or stay constant, the 
evolution of q(t) need not consist of monotonic growth 
from its inflationary value of nearly — 1. Indeed, after 
inflation we know that q(t) decreased from the era of 
radiation domination (q = +1) to the era of matter dom- 
ination (q = and considerable evidence exists that 
it has now dropped below zero 0, 0, H(| . There is no 
reason it could not have experienced such oscillations at 
very early times. The standard formalism for analyzing 
density perturbations breaks down in this case, but our 
techniques should be applicable to any number of horizon 
crossings. 

Note that our formalism depends only upon the scale 
factor a(t), not on the matter theory which supports it. 
Results derived in this way are therefore independent of 
the details of particular inflationary models. For exam- 
ple, it does not matter how many scalars, if any, partici- 
pate in driving inflation. 

It should be possible to compute the one loop effec- 
tive action of the massless, minimally coupled scalar as a 
functional of general a(t). Of course the ultraviolet diver- 
gence is universal and was derived long ago 01 > but the 
ultraviolet finite, infrared contributions have never been 
computed and are, in many ways more interesting. Our 
infrared expression l|143fl for the mode function should 
be of great use there. 

Another application is computing the ultraviolet finite 
part of the coincidence limit of the propagator. Since 
linearized gravitons have the same mode functions, this 
would also apply for quantum gravity. Proceeding in this 
way to compute the infrared part of the non-coincident 
propagator, it should be possible to derive the quantum 
gravitationally induced stress tensor in the presence of a 
class of backgrounds guaranteed to include any solution 
which preserves homogeneity and isotropy. This is an 
essential step in checking whether quantum gravitational 
back-reaction can quench A-driven inflation [141 ] . 

Finally, a simple modification of the formalism devel- 
oped here can be used to derive the mode functions for 
general perturbations in gravity plus a single, minimally 
coupled scalar jij]- These mode functions can be used 
to improve the scalar power spectrum as was done for its 
tensor counterpart 
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